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The fragmentation of diatomic molecules under a stochastic force is investigated both classically 
and quantum mechanically, focussing on their dissociation probabilities. It is found that the quan- 
tum system is more robust than the classical one in the limit of a large number of kicks. The 
opposite behavior emerges for a small number of kicks. Quantum and classical dissociation prob- 
abilities do not coincide for any parameter combinations of the force. This can be attributed to a 
scaling property in the classical system which is broken quantum mechanically. 
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I. INTRODUCTION 

The anharmonicity of molecular vibrations makes the 
dissociation of a molecule by irradiation of laser light a 
relatively difficult task Q]. Consequently, high intensity 
is required for dissociation, for instance, / > 10 15 W/cm 2 
for HF and / > 10 14 W/cm 2 for HCl. At such intensities, 
however, the ionization process dominates and masks vi- 
brational excitation and dissociation. Chelkowsky et al. 
suggested that the dissociation threshold of a diatomic 
molecule can be lowered by two orders of magnitude using 
a frequency chirped laser, and hence dissociation without 
ionization should be possible. In a similar spirit, circu- 
larly chirped pulses have been used by Kim et al. Q 
for the dissociation of diatomic molecules. They found 
that the threshold laser intensity is sufficiently reduced, 
to achieve dissociation without ionization. 

Here, we investigate the possibility for dissociation of 
diatomic molecules under a stochastic force, which could 
eventually be chosen such that ionization is minimal. 
A second motivation for our work is the question, if a 
stochastic driving destroys quantum coherence and even- 
tually brings the quantum and classical evolution close to 
each other. 

We model the force as a sequence of pulses (kicks) 
at random times, each kick carrying an independent 
weight 4] . This type of force, similar to white shot noise, 
has been used to model the passage of ions through car- 
bon foils before Its average strength 7 and the aver- 
age number of kicks (Ar,) determine the dynamics of the 
system which is taken as a Morse oscillator |(| with pa- 
rameters corresponding to Hydrogen Flouride (HF) and 
Hydrogen Chloride (HCl) molecules. Classical and quan- 
tum evolution of the system are conveniently compared 
by using the Wigner transform of the initial wavefunc- 
tion as initial phase space distribution in the classical 
evolution. 

We begin the main text of this paper in section II with 



a brief description of the stochastic Hamiltonian. In sec- 
tion III we explain the classical and quantum method 
with which we solve the stochastic dynamics, which is the 
Langevin equation with test particle discretization and 
Monte Carlo sampling in the classical case and the di- 
rect solution of the stochastic Schrodinger equation with 
a standard FFT Split Operator method with absorbing 
boundary conditions in the quantum case. Results, par- 
ticularly for the dissociation probability will be presented 
and discussed in section IV, while section V concludes the 
paper. 



II. DESCRIPTION OF THE MODEL 

The one-dimensional stochastic Hamiltonian of our 
system reads (atomic units are used unless stated oth- 
erwise) 



H(t) = H - xF(t) = 2- + V{x) - xF(t), (1) 

where the molecular dipole gradient |2| has been ab- 
sorbed into the stochastic force F(t). The Hamiltonian 
Hq describes vibrational motion of the molecule in the 
Morse potential Q 



V(x) = -D e +D e (l- exp(-/3a;)) : 



(2) 



with well depth D e and and length scale (3. The eigenen- 
ergies E n of the Morse oscillator Hq are given by 

E n = Huj e (n+l/2)[l-B(n+l/2)/2], < n < [j] (3) 

where uj e is the harmonic frequency, rib = [j] + 1 is the 
number of bound states with 



J 



1/5-1/2, B = hp(2mD e )~ 1/2 , hoj e = 2BD e . (4) 



The parameters specific to HCl and HF are given in Table 
HI The stochastic force F(t) 0,11 in Eq. (JTJ 
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TABLE I: Parameters of the HF and HC1 molecule for the 
Morse potential. 



stands for a series of random impulses of strength 7 j at 
times ti/\. e F(t) is a kind of white shot noise |7j respon- 
sible for multiple 5— kicks undergone by the molecule, 
where N t is the number of kicks up to time t controlled 
by the Poisson counting process N t . It is characterized 
by the average kicking interval (At) = A -1 about which 
the actual interval Atj = £» — £j— 1 are exponentially 
distributed, similarly as the actual kicking strenghts 7 j 
about their mean 7, 

P(Ati) = \exp(-\Ati), P( 7 <) = 7- 1 exp(- 7 i/ 7 ) . (6) 

We restrict our analysis to positive 7 j and assume that 
7i and ti are mutually uncorrelated random variables gen- 
erated by the distributions functions of Eq. ©. The 
determination of F(t) reduces to the construction of a 
stochastic sequence ( 7 i,ii) which can be done assum- 
ing that the random times ti form a Poisson sequence 
of points leading to a delta correlated process Q. It is 
easy to show 

8] 

that the constructed stochastic force has 
the following properties: 



(F(t)) = 7 A 
(F(t)F(s)) = 2 7 2 A5(i-s)+ 7 2 A 2 

The corresponding power spectrum, i. e. 
transform of (F(t)F(s)), is given by 



7 2 A 



S(u) = 4^= + 7 2 A 2 V2tt5(w). 



(7) 

the Fourier 
(8) 



These properties reveal the difference between the 
present stochastic force (white shot noise) and a pure 
white noise which is delta-correlated with zero mean. 



III. DYNAMICS 
A. Time evolution 

Our system as described is non deterministic due to the 
stochastic nature of its Hamiltonian, but closed. This is 
consistent with a regime of high effective temperature 
and no dissipation. Specifically speaking, the system is a 
simple anharmonic particle which is not coupled to any 
environment (zero dissipation and no diffusion) but sub- 
ject to an external force 0. A perturbative solution for 
this system is in general not possible, because the field 
strengths applied significantly distort the system. We 
are interested in formulation of the dynamics which is 
applicable for the quantum as well as for the classical 



treatment. This can be done in an elegant way by prop- 
agating the Wigner transform W of the density matrix 
t) = \ip(t))(ip(t)\ with a (quantum) master equation 
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Here, L c i and L q represent the classical and quantum 

Liouville operators, respectively, while O stands for the 
superoperator resulting from random kicks undergone by 
the molecule. Unfortunately, solving the master equation 
and constructing O is a complicated task. It is much eas- 
ier to solve the equations of motion derived from Eq. 
directly. 

The classical time evolution obeys the Langevin equa- 
tion 



dp dV 
at ox 



(9) 



while the quantum evolution can be obtained from the 
stochastic Schrodinger equation 



ih 



d \m) 

dt 



H(t)W))- 



(10) 



Both formulations have in common that they must be 
solved over a larger number of realizations N r of the 
stochastic force. Only the average over all realizations 
produces the solution of the Classical Langevin and the 
stochastic Schrodinger equations, respectively. 



B. Initial state 

The molecule is considered to be initially in the ground 
vibrational state |V>o) = |0) with energy Eq (Eq. J3J). For 
the classical propagation we take the Wigner distribution 
of the ground state as initial phase space distribution. 
Analytically, the initial phase space density is given by 

^o( X (o,P) = ^J m e j K_^(o, (id 

where £(x) = 2/B exp(— (3x) and K a is the modified 
Bessel function of the third kind [Tlj . 



C. Classical approach 

The stochastic Langevin equation @ can be solved 
numerically with test particles ("test particles discretiza- 
tion" ) so that the Wigner function is given by 



W(X,p,t) 



Attest 

ACt E 

k=l 



S(x-x k (t))S(p-p k (t)), (12) 



3 



where iVtost is the number of test particles and the 
{xk(t),Pk(t)) are the classically evolved trajectories of the 
test particles. Their initial conditions (xk(0),pk(Q)) are 
Monte Carlo sampled by dividing the phase space into 
small bins Si m 12] . In each Si m the initial conditions 
for ni m test particles are randomly chosen where ni m is 
determined by the value of the Wigner function W(x,p) 
attached to the respective phase space bin, 

= W(x hPm )Si m N test (13) 

with Sim = {xi+l - Xi) X (jp m +l -Pm)- 

For each realization r of the stochastic force Eq. i|12|) 
yields the propagated Wigner function W r {x,p, t) which 
must be averaged over the N r realizations to obtain the 
final result W(x,p,t). 



The dissociation probability for the realization r is then 
given by 

N b -i 

p r d (t) = i- E ( 18 ) 

Classically, P^. c i{t) is given in terms of trajectories 
N + (i) which have positive enery E(t) > at time t. The 
physical result is obtained by averaging Eq. (jTHjl over 
the N r realizations. For the results we will present we 
chose N r = 100 and N tcst — 1000 which was sufficient to 
achieve convergence. 



IV. RESULTS AND DISCUSSIONS 



D. Quantum approach 

For a given realization r, the solution of the stochastic 
Schrodinger equation (|10fl amounts to solve the standard 
time-dependent Schrodinger equation 

IM*)>=EUt,*oM*o)>, (14) 

where U r (t,to) is the evolution operator and \ip(to)) = 
|0). Since the stochastic force consists of instantaneous 
kicks, U r (t,to) can be written as 

U r (t,t )=U {t,t Nr ) Yl exp ( j_xji J U (U+i,U) (15) 

with n r kicks for the realization r and 

U (t i+ i,ti) = exp (-~(t i+ i - ti)Ho^ . (16) 

This representation illustrates how the stochastic driving 
operates. Between two kicks at ti and U+i the molecule 
evolves freely with Uo according to the Hamiltonian Hq 
(Eq. (|16|0 . At each kicking time ti the stochastic force in- 
duces a phase shift by exp (ixjA. In practice, however, 
it is easier to compute \ip r (t)) directly using the Stan- 
dard FFT Split-operator algorithm 13] with absorbing 
boundary conditions. 



E. Dissociation probability 

The observable we are interested in is the quantum dis- 
sociation probability, which is the amount of population 
in the continuum states. However, it is easier to calculate 
the complement, i. e., the population of all bound states 
ip v . It reads for a given realization r 

p:{t) = \mMt))\ 2 - a?) 



A. Quantum and classical dissociation probabilities 

An overview over the results is given in Fig. 1. As a 
general trend one sees that quantum and classical disso- 
ciation probabilities do not coincide, in neither limit of 
the stochastic force (small and large 7 and (At)). Fur- 
thermore, for all parameter combinations, the classical 
dissociation takes off later but increases faster eventu- 
ally overtaking the quantum dissociation and going into 
saturation earlier than in the quantum case. The more 
abrupt classical dissociation can be parameterized with 

P d . tCl ~ - tanh(at + b) + i (19) 

which fits quite well the classical dissociation. The fact, 
that the discrepancy between the quantum and the clas- 
sical curves prevails even in the parameter limits for the 
stochastic force, can be partially attributed to a scal- 
ing invariance. This invariance with respect to the ra- 
tio j /(At) = (F(t)) is obeyed by the classical dynamics 
but not by the quantum dynamics. The scaling invari- 
ance means that for equal average stochastic force (F(t)) 
(compare Eq. (JJJ), the classical dynamics is the same, 
yet on different effective time scales. This can be seen by 
transforming the dynamics to a dimensionless time vari- 
able t = t/(At). The effective Hamiltonian in the new 
time variable r, (H — E)dt/dr, remains invariant against 
changes of (F{t)) = j/(At). 

While the classical dynamics gives qualitatively the 
same picture as the quantum dynamics it does not ap- 
proach the quantum result, not even in the limit of a large 
number of kicks. This is different from a Coulomb system 
under stochastic driving Since it becomes classical 
close to E = (which corresponds formally to the disso- 
ciation threshold here) one can show that the Coulomb 
system itself behaves almost classical, and therefore the 
classical scaling under the stochastic force also applies 
to the quantum Coulomb system close to E = 0. The 
molecular system behaves non-classically, even close to 
the dissociation threshold, which prevents to approach 
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FIG. 1: Quantum and classical dissociation probabilities of HF for 7 and (At) as labelled on each panel, where T or b = 0.085. 
Squares denote the classical dissociation probabilities fit with I I19I . The red (grey) line is a free-particle model, see text. Results 
obtained for HC1 are qualitatively similar. 



the classical scaling under the stochastic force. Interest- 
ingly, the nature of stochastic driving, namely the cancel- 
lation of interferences, does not help to approach classical 
behavior. The reason is that the dynamics in the Morse 
potential without stochastic force differs classically and 
quantum mechanically, particularly for higher energies, 
where the non linearity of the potential is strong. Con- 
sequently, one may ask if under a very strong stochastic 
force, i. e. , without a potential, classical and quantum 
dynamics become similar. 



B. The free particle limit under stochastic driving 

If V = in Eq. (JIJ, i. c. , Ho = p 2 /2m, one sees im- 
mediately with the help of Ehrcnfest's theorem that clas- 
sical and quantum observables should agree since Eq. l[T|l 
only contains linear and quadratic operators. Therefore, 
the state \ip(t)) time evolved under the stochastic driving 
from an initial momentum |^>(Q)) = \p) is simply given 
by \tp(t)) = \p(t)) where p(t) is defined by the classical 
time evolution, starting from the initial momentum pi at 
time tg, 

p t = A 7 (f - t ) + Pi . (20) 

We can define a formal analog to the dissociation proba- 
bility, namely the probability to find a particle after time 



t with positive momentum 

P d (t) = / dp t W( Pl (p t )) , (21) 
Jo 

where Pi{pt) can be obtained from Eq. I|2(J|) and W(p) 
is the initial momentum distribution which we assume 
for simplicity to be Gaussian, 

W(pi) = a/V^cxp(~a 2 p 2 ) . (22) 

Inserting Eq. lt2*21) and Eq. (|2U|) into Eq. lj2*T|) leads to an 
incomplete Gaussian integral with the analytical solution 

P d (t) = |Erfc[-aA 7 (t-t )]. (23) 

We may use this analytical expression (t) with the two 
parameters a, to to fit and interpret the dissociation prob- 
abilities in Fig. 1. At some time to > after a number 
of kicks the systems will have a distribution (with width 
a) of energies and the mean energy may be considered to 
be high enough to approximate the dynamics with a free 
particle Hamiltonian under stochastic driving without a 
potential. As one can see in Fig. 1, this approximation 
becomes increasingly better for larger time in comparison 
to the classical response, while the quantum response re- 
mains different. In Fig. 1-d Pd{t) is plotted for 7 = 0.021, 
A = (AT)" 1 = 235.301, m = 1782.83, a = 0.02135 and 
t Q = 150. 
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V. CONCLUSIONS 

We have proposed and discussed the possibility of dis- 
sociating diatomic molecules by a stochastic force. This 
problem has been explored as function of the characteris- 
tic parameters of the stochastic force, namely the average 
strength of kicks 7 and the average time between kicks 
(At). In view of the effectivity of the stochastic force 
to dissociate the molecule with typical (At) much longer 
than electronic time scales we expect the stochastic force 
to be an efficient way to dissociate a molecule. In con- 
trast to Coulomb dynamics there is no parameter limit 
of the stochastic force where classical and quantum re- 
sults coincide. The reason is the classical scaling of the 
dynamics under the stochastic force which is broken by 
the quantum dynamics. We recall that the present sys- 
tem is a closed one, not coupled to an environment and 
therefore not subject to dissipation and diffusion. For the 
latter case of an open system the classical-quantum cor- 
respondence has been investigated systematically, with 
the general tendency that strong decoherence makes the 



quantum system to behave more classically. In contrast, 
little is known about the quantum-classical correspon- 
dence in the present case of a closed system exposed to 
a stochastic force. 

We hope that our results will stimulate efforts to 
achieve experimental dissociation of diatomic molecules 
using white shot noise. Experiments, using a stochastic 
force similar to the present one, have been successfully 
performed by Raizen and cowerkers [l4| , on the dynamics 
of rotors. 
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